I  A0-A1C7  716 

EFFECTS  OF  VARIABLE  PROPERTIES  IN  FILN  COOLED  TURBULENT 
BOUNDARV  LAVERSIU)  NAVAL  POSTGRADUATE  SCHOOL  HONTEREV 

CA  A  F  HALZ  HAR  86 

1/1 

|  UNCLASSIFIED 

F/S  28/4 

NL 

NAVAL  POSTGRADUATE  SCHOOL 

Monterev,  California 


EFFECTS  OF  VARIABLE  PROPERTIES 
IN  FILM  COOLED 
TURBULENT  BOUNDARY  LAYERS 

by 

Alfred  Frank  Walz  Jr. 

March  1986 

Thesis  Advisor:  P.  M.  Ligrani 


Approved  for  public  release;  distribution  unlimited 


UNCLASSIFIED 


u  st PORT  SECURITY  classification 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  authority 


REPORT  DOCUMENTATION  PAGE 


1ft  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME. OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  ?a  NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School  ^0(je  6g  Naval  Postgraduate  School 


&<  ADDRESS  (City.  Star*,  and  ZIP  Coda) 

Monterey,  California  93943- 5Q00 


7b  ADDRESS  (City.  Star*,  and  ZfPCoda) 

Monterey,  California  93943-5000 


3a  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 


8c  aODRESS  (Cify.  State,  and  ZIP  Code) 


8b  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  epplkebte)  [ 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT  I  TASK 

ELEMENT  NO  I  NO  I  NO 


WORK  UNIT 
ACCESSION  NO 


i '  TITLE  (include  Security  Cnuiticetion) 

EFFECTS  OF  VARIABLE  PROPERTIES  IN  FILM  COOLED  TURBULENT  BOUNDARY  LAYERS 


i.  personal  author(S) 

wait,  Alfred  F.,  Jr. 


’  3a  type  OF  REPORT 

Master's  Thesis 

13b  time  COVERED 

FROM  TO 

is  PAGE  COUNT 

76 

•6  SUPPLEMENTARY  NOTATION 

• 

^  V  >*  i  £r 

' 7  COSATi  CODES  I  iB  SUBJECT  TERMS  (Continue  o~  reverie  it  neceuary  and  identity  by  (Mock  number) 

GROUP  1  $U0GROUP  '^Gas  Turbine  Cooling,  Film  Cooling ,  Turbulent 

Ubundary  Layers ,  Variable  Properties,  C  ' _-f 


'9  ABSTRACT  (Continue  on  reverie  if  neceuary  and  lOentify  by  block  number) 

The  effects  of  variable  properties  on  heat  transfer  in  a  film  cooled 
turbulent  boundary  layer  were  investigated.  A  new  procedure  was  developed 
to  deduce  adiabatic  effectiveness  from  heat  transfer  coefficients  based  on 
the  wall  to  freestream  temperature  difference,  where  both  are  representative 
of  variable  property  flow  conditions.  The  new  technique  was  shown  to  be 
valid  using  data  from  the  literature  for  injection  into  a  turbulent  boundary 
layer  from  one  and  two  rows  of  injection  holes.  From  these  results,  the 
variation  of  the  coolant  to  mainstream  density  ratio  was  shown  to  have  a 
significantly  greater  effect  on  heat  transfer  than  variations  of  viscosity 
and  thermal  conductivity. 


.’0  D  $ TRiSUTlON  /  AVAILABILITY  OF  ABSTRACT  Z1  ABSTRACT  SECURITY  CLASSIFICATION 

QD  lnclassifieo/unlimited  □  same  as  rpt  Q one  users  I  UNCLASSIFIED 


22b  TELEPHONE  (include  Are*  Coat)  22 c  OFM.E  SYMBOL 

f  4  08 )  646-  3382  Code  69XK 


DO  FORM  1473, 84  mar 


8 J  APR  edition  may  be  used  untn  emausted 
All  otner  editions  are  obsolete 


SECURITY  CLASSIFICATION  QC  this  PAGE 

UNCLASSIFIED 


Approved  for  public  release;  distribution  is  unlimited, 


Effects  of  Variable  Properties 
in  Film  Cooled 
Turbulent  Boundary  Layers 


by 


Alfred  F.  Walz  Jr. 

Lieutenant  Commander,  United  States  Navy 
B. S.  ,  Maine  Maritime  Academy,  1977 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  1986 


Author: 
Approved  by: 


7  John  N.  Dyer, 

Dean  of  Science  and  Engineering 


2 


S  V-' •-V-V.  V 


•  •.  v.  .V  >  . 


o'- 


ABSTRACT 
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The  effects  of  variable  properties  on  heat  transfer  in  a 
film-cooled  turbulent  boundary  layer  were  investigated.  A 
new  procedure  was  developed  to  deduce  adiabatic  effective¬ 
ness  from  heat  transfer  coefficients  based  on  the  wall  to 
freestream  temperature  difference,  where  both  are  represen¬ 
tative  of  variable  property  flow  conditions.  The  new  tech¬ 
nique  was  shown  to  be  valid  using  data  from  the  literature 
for  injection  into  a  turbulent  boundary  layer  from  one  and 
two  rows  of  injection  holes.  From  these  results,  the  varia¬ 
tion  of  the  coolant  to  mainstream  density  ratio  was  shown  to 
have  a  significantly  greater  effect  on  heat  transfer  than 
variations  of  viscosity  and  thermal  conductivity.  tTeuto/  ols; 
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NOMENCLATURE 


d  -  Injection  Hole  Diameter 

H  -  Enthalpy 

h  -  Heat  Transfer  Coefficient  with  Film  Cooling 

(spanwise  averaged),  q'V  (TQ0B  -  Tw)  and  q'7  (Treo  -  Tw) 

h£  -  Heat  Transfer  Coefficient  with  Film  Cooling 
(spanwise  averaged),  q"/ (TflW  -  Tw) 

hQ  -  Heat  Transfer  Coefficient  without  Film  Cooling 

(spanwise  averaged),  q'7(T  -  Tw)  and  q’7(Trw  -  Ty) 

2  2 

I  -  Injection  Momentum  Flux  Ratio,  PCU  /p^U^ 

k  -  Thermal  Conductivity 

m  -  Mole  Fractions  of  Constitutant  Gases 

M  -  Blowing  Ratio,  PcUc/paBUao 

P  -  Static  Pressure 

Pr  -  Molecular  Prandtl  Number 

q"  -  Wall  Heat  Flux 

Re  -  Reynolds  Number,  PcTJcd/uc 

T  -  Absolute  Temperature 

U  -  Mean  Velocity 

u  -  Velocity  Component  in  X  Direction 

v  -  Velocity  Component  in  Y  Direction 

W  -  Molecular  Weights  of  Constitutant  Gases 

Distance  from  Downstream  Edge  of  Injection  Holes 


x 


Nondimens ional  Distance 
Thermal  Diffusivity 
Density 

Adiabatic  Film  Cooling  Effectiveness  (spanwise 
averaged)*  eqn.  15  and  22 

Absolute  Viscosity 

Kinematic  Viscosity 

Nondimens ional  Coolant  Temperature,  eqn.  16 

Injection  Angle  Measured  from  Mainstream 
Flow  Direction 

Nondimens ional  Distance  Parameter,  eqn.  25 
Second  Coefficient  of  Viscosity 

rB  SCRIPTS 
Adiabatic  Wall 
Coolant 

Some  Reference  Value 
Total  Condition 
Recovery  Condition 


Freestream 
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Over  the  past  thirty  years,  gas  turbines  have  become 
smaller  and  lighter  with  improved  fuel  economy.  These 
desirable  characteristics  have  been  achieved  by  raising  the 
turbine  entry  temperatures.  To  prolong  useful  blade  life  in 
this  high  temperature  enviornment,  cooling  systems,  thermal 
barriers,  and/or  different  materials  capable  of  withstanding 
higher  temperatures  are  utilizied. 

This  thesis  concentrates  on  film  cooling.  Film  cooling 
protects  components  from  high  gas  temperatures  by  blanketing 
the  component  with  a  film  of  cooler  air.  This  film  insu¬ 
lates  the  component  and  acts  as  a  heat  sink.  The  problem 
encountered  with  film  cooling  is  the  prediction  of  the 
component’s  wall  temperature.  Two  different  convective  heat 
transfer  coefficients  are  currently  used  for  this 
prediction: 

0"  •  VT aw  •  V 


q”  -  hfT  -  T  ) 
n  0™  w' 


Here,  different  temperature  differences  are  used  to  define  h 
and  h£. 

The  objective  of  this  thesis  is  to  find  a  relation 
between  the  two  heat  transfer  coefficients  when  the  effects 
of  variable  properties  are  important. 


[#] 


A.  TURBINE  BLADE  COOLING  FUNDAMENTALS 

For  a  given  size  gas  turbine,  raising  the  turbine  entry 
temperature  results  in  a  higher  specific  power.  This 
enables  the  designer  to  operate  the  ’engine  cycle’  at  a 
higher  pressure  ratio,  resulting  in  improved  thermal  effi¬ 
ciency  of  the  thermodynamic  cycle  and  improved  overall  effi¬ 
ciency.  The  net  result  is  a  lower  specific  fuel 
consumption. 

Raising  the  turbine  entry  temperature  results  in  higher 
component  temperatures.  These  components  become  more  prone 
to  failure  at  higher  temperatures  through  greater  oxidation 
rates  of  their  surfaces,  excessive  creep,  and/or  by  thermal 
fatigue.  Figure  1  (Hennecke,  1982)  shows  a  comparsion  of 
the  rising  turbine  entry  temperature  to  the  rising  blade 
material  temperature  limitations.  Notice  that  uncooled 
metal  components  are  limited  to  a  turbine  entry  temperature 
of  approximately  1200  K.  Current  fuel  economy  goals  require 
a  turbine  entry  temperature  of  1600  K  and  greater. 
Numereous  cooling  systems  have  been  devised  to  maintain 
internal  components  at  or  below  their  maximum  allowable 
temperature.  Vanes,  blades,  disks,  liners,  and  casings  may 
all  require  some  form  of  cooling  protection. 

Figure  2  (Daniels  and  Schultz,  1982)  illustrates  two 
possible  cooling  medium  schemes.  Liquid  systems  use 
external  sprays  and/or  internal  convection  methods.  Air 
schemes  use  external  films  and/or  internal  convection 
methods.  Liquid  systems  have  several  thermodynamic  advan¬ 
tages  over  air  systems.  Dispite  these  advantages,  however, 
air  systems  have  been  universally  accepted  due  to  their 
relative  mechanical  simplicity  to  implement.  Three  prin¬ 
ciple  methods  of  air  cooling  currently  exist.  They  are 


internal  convection,  film,  and  transpiration/effusion. 
Figure  3  (Hennecke,  1982)  illustrates  these  characteristic 
types  of  cooling  for  turbine  blades. 

Internal  convection  is  cooling  by  convection  between  the 
internal  sections  of  the  blade  and  the  air.  This  type  of 
cooling  may  be  subdivided  into  several  methods  consisting 
of: 

( 1)  Straight  through  radial  holes  or  channels  of  a  single 
or  multipass  nature; 

(2)  Enhanced  heat  transfer  due  to  an  enlargement  of  the 
internal  surface  area,  promotion  of  turbulance,  or 
impingement ; 

(3)  An  insert  that  distributes  the  air  according  to  the 
localized  cooling  requirements. 

Film  cooling  is  localized  ejection  of  air  through 
discrete  slots  or  holes  on  the  blade.  This  results  in  a 
film  of  air  around  the  blade  which  is  cooler  than  combustion 
gas  temperature,  acting  as  a  heat  sink.  The  film  also  acts 
like  a  thermal  barrier  between  the  hot  gases  and  the  blade's 
surface.  This  type  of  cooling  may  be  subdivided  into 
several  methods  consisting  of: 

(1)  Single  films;  where  specific  parts  of  the  blade  are 
covered  by  a  film; 

(2)  Full  coverage;  where  many  films  are  placed  close 
enough  together  so  they  cover  part  or  all  of  the 
blade. 

Transpiration  or  effusion  cooling  is  the  bleeding  of  air 
through  the  blade  at  every  point  to  form  a  protective  film. 
The  blade  is  made  of  a  porous  material  to  accomplish  this 
method  of  cooling. 

Typically,  cooling  air  is  extracted  from  various  stages 
in  the  compressor  and  supplied  to  the  components  that 
require  cooling.  Figure  4  (Hennecke,  1982)  illustrates  the 
coolant  flow  path  for  a  single  stage  of  a  typical  gas 
turbine.  Typical  coolant  flow  rates  as  a  percentage  of  the 
gas  turbine’s  total  flow  rate  are  included.  Notice  that  the 
guide  vanes  and  turbine  blades  have  the  highest  coolant  flow 
rates.  These  two  components  are  of  primary  concern  for  they 


have  the  largest  amount  of  surface  area  exposed  to  the  main¬ 
stream  flow  and  consequently  the  highest  thermal  loading. 

Current  applications  use  a  combination  of  the  internal 
convection  and  film  cooling  methods.  Figure  5  (Daniels  and 
Schultz,  1982)  illustrates  a  typical  guide  vane  cooling 
configuration.  Cooling  air  comes  in  from  the  top  and  bottom 
of  the  vane.  Notice  the  internal  convection  methods  used: 

(1)  Air  is  distributed  to  the  internal  surface; 

(2)  Impingement  on  the  surface  after  distribution; 

(3)  Ribs  and  pedestals  used  to.  enlarge  the  surface  area. 
After  the  air  has  cooled  the  vane  by  internal  convection,  it 
is  ejected  into  the  mainstream  and  used  as  a  film  on  the 
exterior  surface.  Initial  stages  utilize  the  internal 
convection  and  film  cooling  principles  as  illustrated  in 
Figure  5.  Latter  stages  may  employ  internal  convection 
cooling  solely  due  to  the  lower  gas  temperatures  encountered 
in  these  stages. 

The  goal  in  using  any  cooling  system  is  to  maintain  the 
internal  components  temperatures  at  or  below  their  maximum 
allowable  temperature.  Achievement  of  this  goal  will  result 
in  high  reliability  and  long  service  life.  In  addition,  it 
is  desirable  to  obtain  (as  closely  as  possible)  a  uniform 
blade  temperature  distribution  to  minimize  thermal  stresses. 

In  the  selection  of  a  cooling  system  to  be  employed,  two 
key  limitations  must  be  kept  in  mind.  These  limitations 
are: 

(1)  The  minimal  amount  of  cooling  air  should  be  used  to 
minimize  it  s  detrimental  effect  on  engine  perform¬ 
ance  ; 

(2)  The  effectiveness  of  a  particular  cooling  method  to 
keep  the  component  sufiiciently  cooled  for  the 
desired  turbine  entry  temperature. 

The  detrimental  effect  on  engine  performance  is  two 
fold.  First,  extracting  air  from  the  compressor  reduces  the 
available  output  power.  Some  work  can  be  recovered  by  rein¬ 
troducing  the  cooling  air  back  into  the  mainstream  flow 


after  cooling.  Second,  ejection  of  air  through  the  blade 
has  a  detrimential  effect  on  the  aerodynamic  efficiency  of 
the  turbine.  Diminished  aerodynamic  efficiency  results  from 
pressure  losses  and  flow  separation  in  the  wakes  trailing 
the  blades. 

Figure  6  (Hennecke,  1982)  shows  the  effectiveness  of 
different  methods  used  to  cool  components  as  a  function  of 
the  turbine  entry  temperature.  ’Cooling  effectiveness’  is 
defined  as  the  ratio  of  actual  to  theoretical  blade  tempera¬ 
ture  reduction.  This  'cooling  effectiveness'  is  not  only  a 
function  of  the  geometric  configuration  involved,  but  also 
the  coolant  flow  rate.  As  the  turbine  entry  temperature 
increases,  so  must  the  required  cooling  effectiveness. 

Figure  7  shows  the  cooling  flow  required  by  various 
cooling  methods  as  a  function  of  the  turbine  entry  tempera¬ 
ture.  Notice  that  convection  cooling  requires  the  most 
cooling  flow  rate  while  effusion  requires  the  least  at  any 
given  turbine  entry  temperature.  Although  this  figure  is 
for  specific  leading  edge  and  coolant  temperatures,  the 
trend  is  the  same  for  other  temperature'  combinations. 

Figures  6  and  7  indicate  that  to  achieve  the  highest 
possible  turbine  entry  temperature  with  minimum  detrimential 
effects,  effusion  cooling  should  be  used.  Effusion  cooling, 
however,  is  not  currently  used  due  to  structural  and  manu¬ 
facturing  difficulties,  and  because  the  porous  blades  often 
become  clogged.  This  is  one  reason  why  film  cooling  methods 
are  currently  employed. 

Logically,  one  could  assume  that  the  more  coolant  pumped 
into  the  film  and  the  fuller  the  coverage,  the  higher  the 
turbine  entry  temperature.  However,  there  are  limitations: 

(1)  The  detrimential  aerodynamic  effects  of  film  cooling 
may  offset  the  gains  of  a  higher  turbine  entry 
temperature ; 

(2)  ’’Lift-Off"  may  occur  where  the  film  ceases  to  protect 
the  surface. 


Thus,  thermal  design  of  vanes  and  blades  with  film  cooling 
must  be  optimized. 

B.  FILM  COOLING 

Film  cooling  is  the  ejection  of  coolant  into  the 
boundary  layer  on  the  blade.  Figure  8  shows  a  typical  film 
cooling  geometry  with  temperature  profiles  for  a  row  of 
holes.  Notice  on  the  top  view,  at  some  distance  downstream, 
the  flows  from  adjacent  jets  converge  to  form  a  continuous 
film  over  the  blade's  surface.  In  this  region,  lateral 
variations  become  negligible,  allowing  a  two  dimensional 
analysis  to  be  conducted  with  reasonable  accuracy. 

The  typical  temperature  profiles  on  Figure  8  illustrate 
the  effect  of  film  cooling  on  heat  transfer.  The  left 
profile  is  without  film  cooling.  Heat  transfer  is  from  the 
freestream  to  the  blade's  surface.  The  driving  force  is  a 
function  of  the  freestream  recovery  temperature  and  wall 
temperature.  The  right  profile  is  with  film  cooling.  Heat 
transfer  is  from  the  freestream  to  the  coolant  and  from  the 
coolant  to  the  blade's  surface.  The  driving  force  is  a 
function  of  the  freestream  recovery  temperature,  coolant 
recovery  temperature,  and  wall  temperature.  As  the  flow 
moves  downstream  from  the  injection  site,  the  coolant  mixes 
with  the  existing  boundary  layer.  Eventually  a  temperature 
profile  similar  to  the  left  will  prevail.  When  this 
happens,  the  heat  sink  and  insulating  effects  of  the  film 
become  minimal. 

Figure  9  illustrates  a  typical  heat  transfer  coefficient 
and  adiabatic  wall  temperature  distribution  around  a  turbine 
blade.  Notice  that  the  highest  temperature  and  heat 
transfer  coefficient  occurs  at  the  nose  of  the  blade,  near 
the  stagnation  point,  and  the  location  of  the  highest 
thermal  loading.  Because  of  this,  injection  sites  are 
located  in  the  nose  and  utilize  the  lowest  coolant  tempera¬ 
ture.  In  the  region  around  the  stagnation  point,  the  heat 


transfer  coefficients  in  Figure  9  are  representative  of  a 
laminar  boundary  layer.  At  approximately  one- fourth  the 
distance  along  the  suction  surface ,  transition  to  a  turbu¬ 
lent  boundary  layer  begins ,  resulting  in  an  increase  in  the 
adiabatic  wall  temperature  and  heat  transfer  coefficient 
approximately  three- fourths  the  distance  along  the  suction 
surface.  The  heat  transfer  coefficients  along  the  pressure 
side  in  Figure  9  are  representative  of  a  turbulent  boundary 
layer  developing  with  a  constant  streamwise  pressure 
gradient  (after  an  initial  drop  from  the  stagnation  point). 
Near  the  mid- chord  position,  the  heat  transfer  coefficients 
begin  to  rise  as  a  strong  favorable  pressure  gradient 
develops. 

C.  BOUNDARY  LAYER  EQUATIONS 

To  predict  the  heat  flux  from  the  mainstream  to  the 
blade,  the  equations  governing  the  physical  situation  must 
be  solved.  Figure  10  illustrates  the  coordinate  system 
used.  Taking  a  differential  element  and  applying  the 

conservation  principles  of  mass,  momentum,  and  energy;  the 
governing  equations  are  derived.  Using  the  governing  equa¬ 
tions,  a  solution  to  the  heat  transfer  problem  is  consider¬ 
ably  difficult  since  flow  and  temperature  distributions  are 
very  complex.  These  equations  are  signif iciently  simplified 
by  making  the  following  assumptions: 

(1)  u  >>  v 

3P 

(2)  —  *  0 

3T  3T 

(3)  —  >>>  — 

3y  3x 

(4)  X  .  .  0 

3 

(5)  No  internal  heat  sources  exist 

(6)  No  chemical  reactions  occur  within  the  element 


(7)  Body  forces  are  negigible 

(8)  No  work  is  performed  by  external  fields 

Applying  these  assumptions  reduces  the  governing  equa¬ 
tions  to  the  boundary  layer  equations.  For  turbulent 
boundary  layers,  the  following  form  is  for  time-averaged 
velocities  and  variable  properties: 
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The  boundary  layer  equations,  as  presented,  are  highly 
nonlinear.  The  momentum  and  energy  equations  are  coupled 
through  the  temperature  dependency  of  density  and  viscosity. 
Thermal  conductivity  is  also  temperature  dependent.  The 


dependence  of  these  transport  properties  on  temperature 
results  in  different  velocity/ temperature  profiles  and 
different  heat  transfer  coefficients  from  those  obtained  by 
assuming  properties  which  are  constant  through  the  boundary 
layer. 

Assuming  constant  properties  exist  through  the  boundary 
layer,  the  boundary  layer  equations  may  be  rewritten  as: 
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The  momentum  and  energy  equations  ((9)  and  (10))  are  uncou¬ 
pled.  This  means  that  the  momentum  equation  can  be  solved 
independently  of  the  energy  equation.  In  addition,  the 
energy  equation  is  linear.  These  two  items  simplify  the 
process  of  solving  the  boundary  layer  equations. 

A  number  of  existing  methods  to  account  for  variable 
properties  in  boundary  layers  without  film  cooling  will  now 
be  discussed. 

Mainstream  gases  consist  of  air  and  the  combustion 
process  by-products  (principally  Carbon  Dioxide,  Carbon 
Monoxide,  and  water  vapor).  Viscosity  and  thermal  conduc- 

tivity  of  this  mixture  can  be  determined  thouqh.  the  semiem— 
pirical  formula  from  White  (1974)  : 
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where 
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This  formula  is  for  a  mixture  of  n  gases.  are  the  molec* 
ular  weights  and  mi  are  the  mole  fractions.  The  exact  same 
formula  is  recommended  for  thermal  conductivity  (y  is 
replaced  by  k).  Figure  11  (White  (1974))  illustrates  the 
temperature  variations  of  viscosity  and  thermal  conductivity 
for  these  gases  at  low  pressures.  Notice  that  for  a  200  R 
temperature  variation,  the  variations  in  viscosity  and 
thermal  conductivity  are  not  very  dramatic.  A  power  law  of 
the  form  (Kays  and  Crawford  (1980)): 


#  -p  \  0 .68 

\  T  / 


(13) 


can  be  used  to  determine  their  temperature  variations.  For 
small  temperature  differences  through  boundary  layers,  say 
20  C  or  less,  viscosity  and  thermal  conductivity  may  each  be 
considered  to  be  constant. 

By  considering  the  mainstream  gases  to  be  a  perfect  gas 
and  from  uniform  static  pressure  through  the  boundary  layer, 
the  density  variation  becomes  inversely  proportional  to  the 
absolute  temperature. 


Comparsion  of  equations  13  and  14  shows  that  temperature 
variations  result  in  larger  density  ratio  changes  than  y  or 
k  changes. 


The  objectives  in  the  thermal  design  of  gas  turbine 
blades  are  to  obtain  an  acceptable  stress  level  and  design 
life.  To  do  this,  a  detailed  blade  temperature  distribution 
must  be  obtained. 

The  thermal  design  of  a  cooled  turbine  blade  is  a 
complicated  energy  balance.  The  design  process  can  be 
broken  down  into  three  major  categories  outlined  below  (in 
sequential  order): 

(1)  Predicting  the  heat  flux  from  the  hot  mainstream 
gases  to  the  blade  s  surface; 

(2)  Performing  a  conduction  analysis  for  the  detailed 
mapping  of  the  blade  s  temperature  distribution; 

(3)  Predicting  the  heat  flux  from  the  blade  to  the 
internal  cooling  air. 

This  thesis  only  addresses  the  first  category:  a  relation 

is  developed  and  verified  between  the  two  types  of  heat 
transfer  coefficient  used  to  predict  wall  temperature. 

A.  CONSTANT  PROPERTY  CONSIDERATIONS 

Equations  1  and  2  are  used  to  calculate  the  blade's  wall 
temperature  after  the  heat  flux  is  obtained  by  solving  the 
boundary  layer  equations. 

1"  "  ht(Taw  ‘  V  (D 

and 
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Both  equations  reflect  the  presence  of  film  cooling. 

In  developing  thermal  boundary  layers  with  film  cooling, 
an  adiabatic  wall  boundary  condition  can  be  used  to  describe 
the  convection  process.  Using  this  concept  reduces  the 
driving  potential  from  a  three  temperature  potential  problem 


to  a  two  temperature  potential  problem.  Equation  1  repre¬ 
sents  this  concept.  Here,  is  an  isoenergtic  heat 
transfer  coefficient.  The  adiabatic  wall  temperature,  Taw, 
is  the  wall  temperature  when  there  is  no  heat  transfer  to 
the  wall.  Adiabatic  wall  temperature  is  usually  expressed 
nondimens ionally  in  terms  of  the  adiabatic  film  cooling 
effectiveness,  n: 


CIS) 


Since  the  entire  surface  is  considered  adiabatic,  T  is  a 

aw 

global  parameter.  As  a  result,  equation  1  is  dependent  on 
the  entire  convection  process  and  the  thermal  boundary 
layer. 

Compared  to  the  large  temperature  differences  between 
the  mainstream  gases  and  the  blade's  surface,  the  surface 
temperature  distribution  can  be  considered  relatively  isoth¬ 
ermal.  An  isothermal  wall  is  ideal  since  thermal  stresses 
are  minimal.  Equation  2  represents  this  concept.  Here,  h 
is  defined  by  the  freestream  stagnation  temperature  and  the 
wall  temperature  difference  for  a  specific  coolant  tempera¬ 
ture.  The  coolant  temperature  is  usually  expressed  nondi- 
mensionally  as  0: 


T  -  T 
rc  '  r® 


There  has  been  considerable  controversy  over  these  two 
methods  of  predicting  the  wall's  temperature.  Two  key 
issues  of  this  controversy  deal  with  how  the  temperature 
difference  is  defined  and  the  effects  of  changing  coolant 
parameters  (temperature,  density,  mass  flow  rate,  etc.).  In 
equation  1,  when  the  temperature  difference  equals  zero, 
there  is  no  heat  transfer.  This  is  not  necessarily  true  if 
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equation  2  is  used.  To  account  for  viscous  dissipation,  the 
freestream  stagnation  temperature,  T0ao,  is  replaced  by  the 
freestreasn  recovery  temperature,  Tra#.  The  heat  transfer 
coefficient,  h,  is  then  redefined  as: 

q”  *  h (T  *  T  ) 

Changing  the  film  cooling  parameters  will  have  different 
effects  on  these  equations.  Generally,  the  recovery  temper¬ 
ature  is  determined  from  freestream  conditions  (temperature, 
velocity,  density),  however,  if  it  is  based  on  near  wall 
conditions,  the  recovery  temperature  in  a  film  cooled 
boundary  layer  is  the  same  as  the  adiabatic  wall  tempera¬ 
ture.  In  equation  1,  both  h£  and  Taw  are  affected,  whereas, 
in  equation  17  only  h  is  affected. 

For  a  given  location,  the  heat  flux  will  be  the  same  for 
either  method.  As  such,  equation  1  is  set  equal  to  equation 
17  and  a  relation  between  the  heat  transfer  coefficients 
exists.  Using  the  definitions  of  q  and  0,  and  dividing  by 
the  heat  transfer  coefficient  for  the  case  of  no  film 
cooling,  then  yields: 


This  equation  relates  adiabatic  and  isothermal  wall  boundary 
conditions  for  properties  that  are  constant  through  the  film 
cooled  boundary  layer.  This  relationship  is  plotted  on 
Figure  12  (Ville,  Cunat ,  and  Richards,  1978).  Note  the 
linearity  between  h/hg  and  9.  This  linearity  is  due  to  the 
fact  that  the  energy  equation  is  linear  for  constant  proper¬ 
ties.  This  equation  is  valid  for  a  given  downstream  loca¬ 
tion,  where  adiabatic  conditions  are  represented  by  the 
horizontal  intercept,  and  isothermal  conditions  by  the  line. 


On  Figure  12,  the  axis  intercept  points  have  special 
meaning.  The  h/ho  axis  intercept  point  represents  the  heat 
transfer  coefficient  when  the  coolant  is  injected  at  the 
same  temperature  as  the  mainstream  (Trc=Tra>,  0=0).  The  9 
axis  intercept  point  indicates  a  wall  condition  that  is 
locally  adiabatic.  At  this  point,  9  equals  the  inverse  of 
the  effectiveness,  n. 

Figure  12  is  only  good  for  one  downstream  location.  As 
you  move  further  downstream,  the  slope  of  the  curve  becomes 
less  and  may  eventually  become  horziontal  indicating  no 
cooling  effect  from  the  film.  This  is  a  result  of  a  rising 
film  temperature  due  to  heat  transfer,  diffusion,  and  mixing 
with  mainstream  gases.  As  the  film  temperature  rises,  the 
value  of  9  decreases  and  the  value  of  n  increases. 


B.  VARIABLE  PROPERTY  CONSIDERATIONS 

With  variable  properties,  the  momentum  and  energy 
boundary  layer  equations  become  coupled  because  of  the 
density  term.  In  addition,  the  energy  equation  becomes 
highly  nonlinear.  Because  of  this,  the  linearity  exhibited 
on  Figure  12  no  longer  exists.  Ligrani  and  Camci  (1985) 
illustrated  this  nonlinearity  that  exists  between  h/h0  and 
9.  They  also  derived  an  empirical  relation  to  relate  the 
constant  property  data  of  Jabbari  and  Goldstein  (1978)  to 
their  variable  property  data. 

A 

To  find  the  relation  between  h/hQ  and  9  for  the  variable 
property  case,  we  start  with  the  defining  equations  for  the 
heat  transfer  coefficients. 

For  variable  properties  through  the  boundary  layer,  the 
variable  property  isoenergetic  heat  transfer  coeffi- 
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cient ,  hp  is  defined  using: 
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The  heat  transfer  coefficient,  h,  is  then  defined  using: 


The  heat  flux  is  now  obtained  by  solving  the  variable  prop¬ 
erty  boundary  layer  equations.  At  a  given  location,  the 
heat  flux  from  equations  19  and  20  will  be  the  same. 
Equating  these  gives: 
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Now  the  adiabatic  wall  temperature  is  expressed  in  terms  of 
an  adiabatic  film  cooling  effectiveness  defined  as: 


(22) 


9  is  defined  as  before  (equation  16).  By  substitution  and 
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rearrangement  the  variable  property  relation  between  h/hQ 
and  6  then  becomes: 


—  =  -41  -  n9) 


(23) 


This  equation  is  for  one  downstream  location  at  one 
value  of  the  density  ratio,  P^P,,,,  and  one  value  of  any  of 
the  following;  blowing  ratio,M;  momentum  flux  ratio,  I;  or 
the  velocity  ratio,  u./u,,,.  A  different  line  results  from 
equation  23  for  each  value  of  the  effectiveness,  n ,  and  the 
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normalized  heat  transfer  coefficient,  h£/hQ.  The  effective¬ 
ness,  n,  may  have  a  strong  dependence  on  the  density  ratio, 
Pc/ p«o»  as  indicated  by  the  single  row  film  cooling  data  of 
Pedersen,  Eckect ,  and  Goldstein  (1977). 

Figure  13  is  a  domain  map  showing  density  ratio  verses 
blowing  ratio  of  the  Pedersen,  Eckect,  and  Goldstein  (1977) 
data.  The  box  symbols  indicate  that  effectiveness  has 
little  or  no  dependence  on  the  density  ratio  for  low  blowing 
ratios.  The  circle  symbols  indicate  that  effectiveness  has 
a  strong  dependence  on  the  density  ratio.  The  triangles 
indicate  that  the  effect  of  the  density  ratio  on  the  effec¬ 
tiveness  and  heat  transfer  is  highly  nonlinear. 


Figure  14  shows  one  possible  effect  of  variable  proper- 
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ties  on  h/hQ  verses  6  data.  Here,  each  of  the  three  linear 
curves  represents  equation  23  with  a  different  n.  The  same 
hf/hQ  was  used  for  each  line  to  simplify  Figure  14.  Now, 
each  line  is  also  for  one  value  of  the  density  ratio,  Pq/Pc, 
which  is  equivalent  to  T^/Tc.  The  variation  of  the  density 
ratio  with  0  is  then  given  by: 
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if  T_ /T_  =  T  /T  for  constant  T  .  The  9  for  a  given 
density  ratio  will  then  be  a  straight  vertical  line 
eminating  from  the  horizontal  axis.  The  intersection  of 
this  vertical  line  and  equation  23  will  be  where  a  given 
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h/1^  data  point  is  located.  The  curved  line  represents  the 
collection  of  the  h/1^  verses  0  data  points.  It  curves 
downward  because  of  the  way  the  data  collects.  As  the 
density  ratio  increases  (vertical  lines  move  to  the  right), 
the  effectiveness  increases.  This  causes  the  0  axis  inter¬ 
cept  point  (1/n)  to  move  to  the  left. 


Using  equation  23,  the  adiabatic  effectiveness  was  used 
to  deduce  h/hQ  using  the  data  of  Pedersen,  Eckect,  and 
Goldstein  (1977)  for  one  row  of  injection  holes.  In  addi- 
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tion,  the  adiabatic  effectiveness  was  deduced  from  h/h  and 
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hf/hQ  from  the  data  of  Ligrani  and  Camci  (1985)  for  two  rows 
of  injection  holes. 

In  the  following  sections,  the  results  are  presented. 

A.  DATA  OF  PEDERSEN,  ECKECT,  AND  GOLDSTEIN  (1977) 

Figure  15  shows  the  Pedersen,  Eckect,  and  Goldstein 
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(1977)  data  used  to  deduce  h/h^  Figure  13  gives  the 
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blowing  ratio  and  the  density  ratio  values  for  which  n  was 
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selected  for  deducing  h/hQ. 

A  value  of  1.  0  was  assumed  for  the  h/hQ  axis  intercept. 
Although  in  many  cases,  a  reasonable  assumption,  Metzger, 
Carper,  and  Swank  (1968)  showed  that  this  axis  intercept  may 
not  be  equal  to  1.0.  They  showed  that  the  injection  has  an 
effect  on  the  heat  transfer  coefficient  even  when  the 
coolant  is  injected  at  the  freestream  temperature,  primarily 
due  to  boundary  layer  thicking. 

Normalized  heat  transfer  coefficients  as  a  function  of  9 
are  shown  on  Figures  16  thru  19  for  the  Pedersen,  Eckect, 
and  Goldstein  (1977)  data.  Four  axes  are  used  to  denote 
four  downstream  locations  (x/d).  The  light  linear  lines 
represent  equation  23  for  a  different  value  of  effective- 
ness,  n.  The  dark  lines  represent  the  collection  of  h/1^ 
verses  9  data  points. 

Figure  16  is  for  a  blowing  ratio,  M,  of  0.213.  This 
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figure  shows  that  h/tig  verses  9  is  steeper  at  small  x/d  and 
becomes  flatter  as  the  distance  from  the  injection  site 
increases  (increasing  x/d).  This  figure  shows  reasonable 
linearity  with  some  nonlinearity  near  the  injection  site. 


Figure  17  is  for  a  blowing  ratio,  M,  of  0.515.  This 
figure  shows  the  same  trends  as  with  Figure  16.  Only  now, 
the  nonlinearity  near  the  injection  site  is  more  pronounced 
and  is  also  evident  at  a  location  further  downstream. 

Figure  18  is  for  a  blowing  ratio,  M,  of  1.05.  This 
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figure  shows  that  the  variation  of  h/hQ  with  0  is  nonlinear. 
As  the  distance  from  the  injection  site  increases 

A 

(increasing  x/d),  the  h/hQ  variation  with  0  becomes  more 
horizontal.  Notice  that,  in  general,  at  a  particular 
blowing  ratio,  the  slope  becomes  more  negative  as  0 

increases. 

Figure  19  is  for  a  blowing  ratio,  M,  of  1.96.  This 

figure  shows  a  highly  nonlinear  behavior  due  to  the  ’lift¬ 
off'  of  film  cooling  jets  from  the  wall,  and  possible  pene¬ 
tration  through  the  boundary  layer. 

B.  DATA  OF  LIGRANI  AND  CAMCI  (1985) 

Figure  20  (Ligrani  and  Camci,  1985)  illustrates  the 
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downward  sloping  trend  for  variable  property  h/hQ  verses  0 
data.  This  figure  also  shows  the  linear  h/hQ  verses  0  data 
(constant  property  case).  Notice  the  similarity  between  the 
variable  property  trend  here  and  the  trends  of  the  Pedersen, 
Eckect,  and  Goldstein  (1977)  data  on  Figure  18. 
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Using  h/hQ  and  hf/h0  from  Ligrani  and  Camci  (  1985),  n 
was  deduced  from  equation  23.  Results  are  plotted  as  a 
function  of  x/d  in  Figures  21  thru  31,  and  as  a  function 
of  £  in  Figures  32  thru  42.  The  later  function,  E,  is  a 

downstream  distance  normalized  by  the  blowing  ratio.  It  is 
defined  by: 
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where  S=itd/6  is  the  width  of  an  equivalent  two  dimensional 
slot  (Jabbari  and  Goldstein,  1978)  .  Goldstein  (1971)  gives 


two  equations  to  predict  laterally  averaged  effective¬ 
ness,  n»  for  two  dimensional  film  cooling.  These 

predictions  are  given  by: 

“  -  (1  ♦  0.2495)'0*8  (26) 

and 


where 

6  -  1  ♦  l.S(10'4) (Rec(^£))sinT  (28) 
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These  predictions  are  plotted  on  the  n  verses  £  graphs 
(Figures  32  thru  42). 

Figures  21  thru  42  are  for  different  blowing  ratios  (M  = 
1.242,  1.003,  0.481,  0.661,  1.557),  momentum  flux  ratios  (I 
=  1.382,  0.925,  0.417),  and  velocity  ratios  (^/u^  *  0.630, 
1.  113,  0.922).  In  general,  these  figures  show  that  as  the 
blowing  ratio,  momentum  flux  ratio,  and  velocity  ratio 
increase,  the  adiabatic  film  cooling  effectiveness 
increases.  The  largest  increase  in  n  occurs  for  the  highest 
value  of  0.  This  corresponds  to  the  case  of  the  highest 
density  ration.  Oiu  Figures  21  thru  31,  the  crossover  point 
between  the  01  and  02  curves  moves  downstream  (from  x/d  =  20 
to  x/d  =40)  as  the  momentum  flux  ratio  and  velocity  ratio 
increase.  On  Figures  32  thru  42,  as  the  blowing  ratio, 
momentum  flux  ratio,  and  the  velocity  ratio  decrease,  the 
data  moves  to  a  higher  value  and  spreads  out.  All  data 
points  on  the  n  verses  £  graphs  are  below  the  predictions  of 
equations  26  and  27  but  do  follow  the  same  trend  as  the 
predictions. 


V.  SUMMARY 


A  relation  was  derived  that  relates  the  normalized  heat 
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transfer  coefficients  h  and  hr.  These  heat  transfer  coeffi- 
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cients  are  for  isothermal  wall  boundary  conditions,  h,  and 
adiabatic  wall  boundary  conditions.  h£.  Each  are  for  when 
variable  properties  are  considered  to  exist  through  the 
boundary  layer.  This  relation  was  presented  as: 


When  the  effects  of  variable  properties  are  important  in 

film  cooled  turbulent  boundary  layers,  the  variable  property 

/\ 

effectiveness,  n,  and  the  above  equation  are  valid  for  one 
density  ratio  and  blowing  ratio  only. 

The  normalized  heat  transfer  coefficient,  h/ho,  was 
deduced  from  the  Pedersen,  Eckect,  and  Goldstein  (1977)  data 
on  injection  with  one  row  of  holes.  The  adiabatic  film 
cooling  effectiveness,  n,  was  deduced  from  the  Ligrani  and 
Camci  (1985)  data  on  injection  with  two  rows  of  holes.  The 
data  deduced  from  equation  23  showed  the  nonlinear  trend  of 
h/h  verses  8  and  as  the  density  ratio,  p  / pB,  varied. 

The  effect  of  density  ratio  on  h/ho  verses  0  results  in 
a  change  of  the  characteristic  curve  from  linear  to  nonli¬ 
near.  This  change  is  not  always  evident.  It's  existance  or 
extent  is  dependent  upon  the  injection  parameters  (blowing 
ratio,  momentum  flux  ratio,  and  velocity  ratio)  and  the 
downstream  distance.  This  effect  is  due  to  the  nonlinearity 
introduced  into  the  energy  equation  and  the  coupliing  of  the 
momentum  and  energy  equations  when  variable  properties 
through  the  boundary  layer  are  considered. 


Experiments  need  to  be  conducted  to  determine  the  accu¬ 
racy  and  dependability  of  the  proposed  variable  property 
relation  on  injection  system  geometry,  injection  flow  rate 
parameters,  fluid  mechanics  of  the  flow,  the  thermal 
boundary  layer  conditions,  and  other  parameters.  In  order 
to  completely  verify  the  proposed  equation,  simultaneous 
measurements  of  variable  property  effectiveness  and  variable 
property  heat  transfer  coefficients  (based  on  the  wall  and 
recovery  freestream  temperature  difference)  need  to  be 
conducted. 
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1.  Internal  Cooling 

-  straight  through 
or  multipass 

-  enhanced  heat  transfer 

-  insert 


2.  Film  Cooling 


-  single  films 

,<F 

-  full  coverage 

Effusion  Cooling 

*  tag : 

Fig.  3  Characteristic  Types  of  Blade 
Cooling  Configurations. 


Fig.  4  Coolant  Flow  in  a  Typical 
High  Pressure  Turbine  Stage. 
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Fig.  5  Coolant  Flow  in  a  Typical 
Nozzle  Guide  Vane. 


Effusion  Cooling 
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Cooling  Effectiveness 


6  Turbine  Entry  Temperature  Verses 
Required  Cooling  Effectiveness. 


Fig.  9  Distribution  of  Heat  Transfer  Coefficient  and  Adaibatic  Wall 

Temperature  Around  a  Typical  Blade. 
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Variation  of  n  with  x/d  for  M=0. 66 
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